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Abstract. We report on near-infrared medium-resolution spectroscopy of a sample of luminous and ultra-luminous 
infrared galaxies (LIRGs-ULIRGs), carried out with SOFI at the ESO 3.5m New Technology Telescope. Because 
of wavelength dependence of the attenuation, the detection of the Paa or Br7 line in the Ks band should provide 
relevant constraints on SFR and the contribution of an AGN. 

We find, however, that the intensities of the Paa and Br7 lines, even corrected for slit losses, are on average 
only 10% and 40%, respectively, of that expected from a normal starburst of similar bolometric luminosity. The 
corresponding star formation rates, after correcting for the attenuation derived from the NIR-optical emission line 
ratios, are 14% and 60% of that expected if the far infrared luminosity were entirely powered by the starburst. 
This confirms the existence of a recombination photon deficit, particularly in the case of the Paa line, already 
found in the Br7 line in other infrared galaxies of similar luminosity. 

In discussing the possible causes of the discrepancy, we find unlikely that it is due to the presence of an AGN, 
though two objects show evidence of broadening of the Paa line and of the presence of coronal line emission. In 
fact, from our own observations and data collected from the literature we argue that the studied galaxies appear 
to be predominantly powered by a nuclear starburst. 

Two scenarios compatible with the present data are that either there exists a highly attenuated nuclear star 
forming region, and/or that a significant fraction (~80% ) of the ionizing photons are absorbed by dust within 
the HII regions. We suggest that observations in the Bra spectral region could constitute a powerful tool to 
disentangle these two possibilities. 
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1. Introduction 

The advent of the Infrared Space Observatory (ISO) com- 
bined with the availability of new ground based facilities 
such as SCUBA on the JCMT have resulted in the discov- 
ery of numerous distant galaxies with enhanced infrared 
(IR) emission (e.g. Elbaz et al. [I999. Small et al. 2000:, 
Barger et al. 120001 Gear et al. I20UUL Franceschini et al. 
IMTTl Small et al. 110021 Ivison et al HUU^ . 
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* based on observations collected at the European Southern 
Observatory, Chile, ESO No. 67.A-0593, 71.A-0707. 



It is not yet clear what mechanisms power the bolo- 
metric luminosity of these luminous (LIRGs) and ultra 
(ULIRGs) luminous infrared galaxies, with bolometric lu- 
minosities exceeding lO^^i© and IO^^Lq, respectively. 

Perhaps, the tightest constraint on the nature of the in- 
frared luminosity is that they mostly obey the FIR/Radio 
correlation of normal star forming galaxies (Sanders & 
Mirabel 1996) . However the presence of broad emission 
line components, a certain degree of polarization, power- 
law near infrared colours, warm far-infrared (FIR) spectra, 
high radio brightness temperature in the milli-arcsecond 
central structure and, finally, hard X-ray emission are of- 



2 



J.R. Valdes et al.: NIR Spectroscopy of LIRGs 



ten detected, suggesting that a significant contribution to 
the luminosity may come from mass accretion onto a com- 
pact central object. 

These galaxies are often found to be interacting sys- 
tems, indicating that the two phenomena may be triggered 
at once by dynamical interaction (e.g. Rigopoulou et al. 

Optical studies have been performed in order to estab- 
lish the nature of these objects (e.g. Veilleux et al. 1999 ), 
but it is difficult to disentangle between AGN and star for- 
mation processes, because the bulk of their luminosity is 
strongly attenuated by dust and the re-emitted FIR spec- 
trum keeps little memory of the primary power source. 

Even in objects without clear signatures of the pres- 
ence of AGN, there are significant discrepancies between 
different observables. For example Poggianti, Bressan & 
Franceschini H2(JU1|I . in analyzing a sample of Infrared 
Galaxies with log(L7j^/L0)>11.5, concluded that the Star 
Formation Rate (SFR) detected from optical diagnostic 
tools (even corrected for attenuation) amounts to about 
1 /3 of that inferred from the FIR. However, if such a large 
fraction of stellar activity is hidden by dust, one may won- 
der whether also the nuclear activity is invisible at optical 
wavelengths. 

A natural way to avoid strong dust obscuration is to 
look in the near infrared where, being the attenuation only 
a fraction (~l/7) of that in the optical, one may hope to 
obtain an unbiased estimate of the SFR and at the same 
time to search for the possible presence of an obscured 
AGN by means of broadening of permitted emission line 
and/or presence of coronal lines. 

With these purposes we begun a medium resolution 
NIR spectroscopic study with SOFI at the ESQ 3.5m New 
Technology Telescope, of a sample of ULIRGs selected 
from Genzel et al. I|1998l) and Rigopoulou et al. (1999) 
samples. In this paper we present and discuss these new 
observations. 

The plan of the paper is as follows. In Sections |21 we 
present the observations and data reduction, with an anal- 
ysis of the emission line properties of target galaxies pre- 
sented in subsections 12 . II and 12 . 21 All the process to calcu- 
late the Star Formation rates are discussed in Section |31 
In particular, the technique we adopted to calculate the 
slit losses and to evaluate the extinction are presented in 
subsections 13.11 and 13.21 In Section 0] we show that, even 
correcting for extinction and losses from the slit, there is 
a large Paa/IR deficit, and we discuss its possible origin. 
Our conclusions are summarized in Sectional We assume 
Ho=65 [km/s/Mpc], 17m=0.3, r2A=0.7 throughout. 

2. Observations and Data Reduction 

The selected galaxies were drawn from the IRAS 1.2Jy 
sample (Fisher et al. 11995(1 . and from Genzel et al. 
H1998|l and Rigopoulou et al. 1(1999(1 . The sample galaxies 
cover the interval of solar luminosities 11.45<log(L/fl/L 
0)<12.73. The IR and radio properties of our target galax- 
ies are listed in TableQ The penultimate column indicates 



the FIR/Radio correlation parameter at 1.4GIIz, (see eq. 
^ below), which for star forming galaxies is known to be 
q~2.35±0.15 over a wide range of infrared luminosities 
(e.g. Sanders & Mirabel [TMH| . 

Medium-resolution Ks-band NIR spectroscopy of the 
galaxies was performed with the SOFI infrared spectrom- 
eter (Moorwood et al. 11998(1 at the Nashmyth A focus of 
the ESO 3.5m New Technology Telescope (NTT), during 
August 2001 and August 2003. During the 2003 run, 5 
minutes of Ks-band imaging were obtained as well. 

The instrumental set consisted of the medium resolu- 
tion grism (pixel scale — 0.292 arcsec), with a dispersion 
of 4.62 A/pixel and a spectral coverage between 2.0 /im 
and 2.3 /im. The Ks filter works as order sorting filter. 
We adopted 1 arcsec slit, providing a spectral resolution 
of R = A A/A « 2000 (Aw w 150 km s^^) at 2 /im. 

Data were acquired by placing the galaxy at two differ- 
ent positions along the slit, through the standard nodding 
technique; the nod throw was set to 60 arcsec. 

In some cases, given the peculiar structure of the 
galaxy, the slit has been oriented with two different posi- 
tion angles (see Table 

The standard reduction pipeline for long-slit NIR 
spectra was performed by using routines of the Image 
Reduction and Analysis Facility (IRAF^). The spectra 
were fiat-field corrected in the usual manner and sky was 
subtracted through the nod on slit procedure. 

Wavelength calibration and slit-curvature correction 
was obtained with the Xenon/ Neon internal lamp and 
checked on sky IR emission lines. 

Corrections for telluric features and flux calibration 
have been performed by observing nearly featureless hot 
B and/or solar- like telluric and spectro- photometric stan- 
dard stars, from the Hipparcos catalogue. Atmospheric 
transmission variations were compensated by observing 
the standard stars contiguously, at similar (as near as pos- 
sible to target galaxies) air masses. The adopted procedure 
consists in dividing the raw galaxies spectra by that of the 
standard star and then multiplying by its intrinsic emit- 
ted spectrum. The latter has been obtained by normalizing 
Pickles ((1998(1 stellar templates of the proper spectral type 
to the V-to-K observed magnitudes of the standard stars. 
In the case of not optimal seeing conditions, we corrected 
the raw spectra of the standard stars for slit-losses, before 
applying flux calibration. 

The exposure time, the position angle (PA) of the slit, 
the spectral type of the atmospheric calibrators and their 
air masses with respect to those of the observed galaxies 
are given in Table El 



^ The package IRAF is distributed by the National Optical 
Astronomy Observatory which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 



J.R. Valdes et al.: NIR Spectroscopy of LIRGs 




(s)s 


1 1 1 1 1 1 1 1 1 1 1 1 
^ IRASO0 188-0856 


1 

w 


X 



2.15 2.2 



IRAS22491-1808 




1.9 1.95 2 2.05 2.1 2.15 





IRAS1 9335-3632 



2.05 2.1 1.85 1.9 1.95 



2.05 2.1 



IRAS23230-6926 



; 0.1 - 



1.9 1.95 



IRAS23389-6139 




1.85 1.9 1.95 



Fig. 1. Spectra of the sample galaxies showing emission in only one nucleus. No slit-losses, nor extinction correction have been 
applied. The identified emission features are highlighted. Fluxes are in units of [10~^^ergs~^cm~^A~^]. 



2.1. Emission line properties 

NIR spectra of galaxies in our sample, presented in Figures 
Hand 121 are characterized by strong Paa or Br7 emission, 
depending on their redshift. In the Paa region we also 
detect weak II2 lines from higher order vibrational transi- 
tions (H2 1-0S(5) at 1.835//m, H2 1-0S(3) at ^.957^lm, and 
II2 1-0S(2) at 2.033/im). Almost all these galaxies show a 
contamination of the Paa blue wing by the Hel emission 



at 1.868/im, and in some of them the weak JiiS emission at 
1.945/im is detected. Galaxies observed in the Br7 region, 
show H2 lines from lower order vibrational transitions (II2 
1-0S(3) at 1.957^m, H 2 1-0S(2) at 2.033 ^^m, H2 l-OS(l) 
at 2.121^m, and H 2 l-OS(O) at 2.223/im). It is to note the 
relatively strong Hel emission at 2.058/im in the case of 
IRAS 23128-5919. 
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Fig. 2. Spectra of the sample galaxies showing emission in two nuclei. No sht-losses, nor extinction correction have been applied. 
The identified emission features are highlighted. Fluxes are in units of [lO^^^erg s^^cm ^^A^^]. 



Table 1. Infrared and radio properties of the observed galaxies 



Object 




z 


D 








SlOOfim 


S25/S60 




Sir 




SRadio 




L/C " 








[Mpc] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 






[W m^] 




[mJy] 




7.7^im 


IRAS00085-1223 





0198 


80.3 


0.395 


2.372 


16.62 


16.97 


0.143 


1 


39x10" 


la 


66.9" 


2.48 




1RAS00188-0856 





1284 


648.6 


0.117 


0.372 


2.59 


3.40 


0.144 


2 


44x10" 


13 


15.7" 


2.36 


1.53 


1RAS00582-0258 





0874 


373.7 


0.110 


0.341 


1.21 


0.99 


0.282 


1 


33x10^ 


13 


10.1" 


2.14 




IRAS01077-1707 





0334 


158.4 


0.301 


0.846 


6.48 


10.4 


0.130 


6 


40x10" 


13 


43.8" 


2.32 




1RAS02411+0354 





1436 


638.1 


0.085 


0.224 


1.37 


1.94 


0.164 


1 


40x10" 


13 


6.6" 


2.44 


3.17 


1RAS06206-6315 





0924 


455.9 


0.069 


0.294 


3.96 


4.58 


0.074 


3 


10x10" 


13 


21.9" 


2.53 


3.69 


IRAS 19335-3632 





0821 


349.8 


0.107 


0.154 


1.19 


1.78 


0.129 


1 


27x10" 


13 


7.1" 


2.36 




IRAS20100-4156 





1295 


655.0 


0.135 


0.343 


5.23 


5.16 


0.065 


4 


00x10" 


13 


20.3^^ 


2.49 


1.92 


1RAS22206-2715 





1314 


579.1 


0.096 


0.160 


1.75 


2.33 


0.091 


1 


61x10" 


13 


6.3" 


2.56 




1RAS22491-1808 





0777 


379.7 


0.119 


0.549 


5.44 


4.45 


0.101 


4 


13x10" 


13 


5.9" 


3.02 


2.85 


1RAS23128-5919 





0446 


184.7 


0.250 


1.590 


10.80 


10.99 


0.147 


9 


07x10" 


13 


37.3" 


2.72 




1RAS23230-6926 





1063 


529.2 


0.058 


0.295 


3.74 


3.42 


0.079 


2 


77x10" 


13 


32.0" 


2.31 


1.50 


1RAS23389-6139 





0927 


457.5 


0.063 


0.244 


3.63 


4.26 


0.067 


2 


83x to- 


13 


163.9" 


1.62 


1.33 



a) 843MHz flux from Mauch et al (2003); b) 1.4GHzflux from NRAO VLA Sky Survey (Condon et al 1998); c) Condon et al 
(1996) 

d) FIR/Radio correlation parameter (Sanders & Mirabel 1996). 843MHz fluxes were converted to 1.4GHz fluxes by assuming 
F{v) oc 

e) L/C from Lutz et aL lB^ 



We find a good correlation between the Paa and IR 
luminosity, as shown in Figure 0] The dashed line in the 
Figure is the zero-intercept linear correlation we obtain, 
log(Lpaa/L/fl) = -4.69. IRAS 02411+0354 has been ex- 
cluded from the analysis because it falls significantly off 
the correlation. In the case of Br7 luminosity, we have 
only three galaxies, but our values fall on top of the re- 
lation found by Goldader et al. ()1997| . The dashed line 



is a zero-intercept linear correlation with slope quoted by 
Goldader et al. (1997), log{h Br ffi^m) = -4.96. 

The spatial distribution of the Paa emission in IRAS 
00582-0258, IRAS 23230-6926 and IRAS 23389-6139 ap- 
pears compact. The peak of the Gaussian profile of the 
line emission corresponds to the peak obtained from the 
continuum emission alone, coming from the nucleus of the 
galaxies. In contrast, IRAS 22491-1808 shows a broader 
spatial distribution in Paa emission, with some minor 
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Table 2. Observational Parameters 



IRAS name 


Integ. 

(see') 


P.A. 

fdee) 


Atmos. 
Calib. 


Gal./Std. 

Air Mass 


00085-1223 


3600 


38 


B9V 


1.11/1.11 


00188-0856 


7200 





F8V 


1.13/1.46 


00582-0258 


3600 


55 


GOV 


1.18/1.17 




3600 


67 


B9V 


1.11/1.16 


01077-1707 


480 


30 


F6V 


1.22/1.25 


02411+0354 


6600 


-54 


GOV 


1.22/1.07 




6600 


90 


GOV 


1.19/1.07 


06206-6315 


2400 


145 


GOV 


1.57/1.69 


19335-3632 


5400 


43 


B9V 


1.19/1.47 


20100-4156 


7500 


25 


G2V 


1.32/1.05 




8400 


-56 


B9V 


1.11/1.12 


22206-2715 


3600 


-32 


GOV 


1.11/1.04 




3600 


61 


GOV 


1.37/1.04 


22491-1808 


7200 


74 


G3V 


1.02/1.08 


23128-5919 


4800 


35 


B4V 


1.32/1.10 


23230-6926 


4500 


158 


F7V 


1.45/1.51 


23389-6139 


480 


8 


F5V 


1.20/1.03 



peaks superposed on a broader single Gaussian profile 
along the slit. 

Six galaxies, IRAS 00188-0856, IRAS 02411-^0354, 
IRAS 19335-3632, IRAS 20100-4156, IRAS 22206-2715, 
and IRAS 22491-1808 show double profiles, corresponding 
to the double nuclei observed in the corresponding NIR 
images. Finally, IRAS 06206-6315 is best fitted with three 
Gaussian profiles. This complexity is associated with the 
bright regions seen in the H and I images (Bushouse et al. 

mm . 

In the case of galaxies observed in the Br7 spectral do- 
main, IRAS00085-1223 and IRAS23128-5919 show a sin- 
gle Gaussian profile, while IRAS 01077-1707 is best fitted 
with three Gaussian profiles. 

The observed emission line properties of target galaxies 
are shown in Table |3| Only for those galaxies with suffi- 
cient S/N have we considered the emission properties of 
multiple components separately while, in all other cases, 
we have folded the separate contribution together. The 
FWHM velocities have been corrected by instrumental re- 
sponse after subtracting in quadrature the instrumental 
FWHM. 

2.2. Detection of AGN signatures 

Typical signatures of central non thermal activity are the 
presence of broad emission line components and/or high 
excitation lines. 

None of the objects observed in the Br7 region show a 
broad component of the Br7 line. In the case of the Paa 
galaxies, we have attempted to isolate a broad compo- 
nent of the Paa emission, by forcing the fit of the shape of 
the line with two Gaussians. With this procedure, in IRAS 
00188-0856 and IRAS 00582-0258 we could detect a broad 
component of FWHM ~ 2339 and 2210 km/s, respectively. 
Figure|3|shows the results of the multi Gaussian fitting for 



these galaxies. In the case of IRAS 00188-0856 we have 
also taken into account the contribution of the Hel line at 
1.8689/im to the general shape of the line. The broad com- 
ponent corresponds to about 42% and 80% of the Paa flux 
in the narrow component for IRAS 00188-0856 and IRAS 
00582-0258, respectively. The Paa intensities reported in 
Table 13 correspond only to the narrow component. 

Another AGN diagnostic in the NIR, is the presence of 
the high excitation line [SiVI]A1.9628/im. This diagnostic 
is particularly useful in Sy2 galaxies, that lack broad-line 
components. The [SiVI]A1.9628/im line falls slightly long 
wards of the H2 1-0S(3) 1.957 //m line, causing an asymme- 
try in medium resolution spectra (e.g. Vanzi et al. I2002|l . 
While our spectral resolution allow us to disentangle these 
two lines, the S/N is not sufficient to draw firm conclu- 
sions. However, [SiVI]A1.9628/j,m seems to be present in 
the spectra of IRAS 00188-0856 and IRAS 00582-0258, the 
same galaxies for which we have hints of a broad emission 
component at Paa. In all other cases there is no evidence 
of this coronal line. 

3. The star formation rates. 

In this section we compare the star formation rate derived 
from the Paa or Br7 emission lines, with that deduced 
from the IR, assuming that both originate in an ongoing 
vigorous starburst. 

In order to perform a proper analysis of the star for- 
mation rate from the emission lines we need to correct the 
measured fluxes for slit losses and for dust extinction. 

3.1. Slit loss corrections 

When comparing different observations, it is important 
to evaluate any discrepancy caused by the different aper- 
tures. If the source has a certain degree of symmetry, 
it is possible to extrapolate the mono dimensional spa- 
tial proflle observed along the slit, to a two-dimensional 
surface brightness distribution, and to estimate the flux 
that would be recevied by a given aperture. We deflne as 
"slit losses" the ratio between the flux corresponding to a 
given aperture and the flux received within our rectangu- 
lar aperture. 

In our galaxies, the two-dimensional flux distribution 
was obtained by rotating around the nucleus the best fit 
Gaussian profile of the spatial distribution of the Paa or 
Br7 line. If the slit covers more than one nucleus, dif- 
ferent Gaussian profiles were fitted for each nucleus. As 
an example. Figure [S] shows the case of IRAS20100-4156, 
at two position angles P.A. =-24 and 56 respectively. Slit 
loss coefficients, Chq and Cm, needed to compare our Paa 
fluxes to Ha and IR fluxes taken from the literature are 
reported in Tabled] 

3.2. Attenuation 

Correction for attenuation can be derived by comparing 
the observed hydrogen line emission intensity ratios, Robs? 
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Fig. 3. Possible detection of a broad emission component in the IRAS 00188-0856 (left panel) and IRAS 00582-0258 
(right panel). The width of the broad component is FWHM~2339 km/s for IRAS 00188-0856, and FWHM~2210 km/s 
for IRAS 00582-0258. 



log{Lp„/L,„) = -4-69 



los(LBr/L„) = -4.97 
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Fig. 4. Relation between Paa and IR luminosity (upper 
panel) and Br7 and IR luminosity (Lower panel). Dashed 
line in the upper panel is our linear relation L(line) = 
Const X L(IR), while that in the lower panel has ben 
taken from the value quoted by Goldader et al. 1997. 



with the intrinsic values predicted by models of nebular 
emission, Ri„t. 

However, since the only strong Hydrogen line present 
in our spectra is either Paa or Br7 we need to compare our 
fluxes with other observations taken from the literature. 
Due, Mirabel & Maza l)1997|l . Veilleux, Kim & Sanders 
H1999|l . and Veilleux et al. H1995|l provide optical spec- 
troscopy within slit apertures of 1.25", 2 kpc and 4 kpc re- 
spectively, for almost all galaxies in our sample. Observed 
intensities of the Ha and H/3 emission lines are reported 
in columns two and three of Table El 



120100-4156" 
PA =-34 



.VA. Va. 




Fig. 5. Example of calculation of slit losses for IRAS 
20100-4156. A Gaussian fit of the Paa spatial profile of 
IRAS 20100-4156 is performed along the slit, for the two 
nuclei when needed. The 1" slit width is reported as 
shaded. The Gaussian is thus rotated around each nucleus 
and the contributions within the slit aperture and a suit- 
able circular aperture around each nucleus, are separately 
evaluated. 

There exists a wide literature concerning predicted 
intrinsic line ratios of Hydrogen recombination lines. In 
general these values depend only very slightly on the as- 
sumed density of the emitting region (usually between 
10^ — lO^cm^"^) but they depend more on the electronic 
temperature. The latter may be assumed or may be the re- 
sult of detailed modeling of HII regions, accounting for the 
hardness of ionizing spectrum, geometry and metallicity 
of the gas. To investigate this point, we have run a set of 
CLOUDY (Ferland, 2003) models for different HII regions, 
characterized by different ages of the ionizing stellar pop- 
ulation, different metallicity and gas density. The intrinsic 
ratio decreases at increasing metallicity (lower electronic 
temperature) and increasing age (lower hardness, neglect- 
ing a hardening of the spectrum due to the presence of 
Wolf Rayet stars). For example, for the ratio R fj^'^""" we 
find values as high as 8.5 at Z=0.008 and young ages and 
as low as 7.0 in metal rich HII regions and/or old clus- 
ter ages. Furthermore, more realistic galaxy models must 
take into account a law of star formation (almost contin- 
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IRAS name 




Line 


At, 


Flux 




EW 


FWHM 


FWHM 








[/iin] 


[org s cm ^] 


[restA] 


[restA] 


[km-i] 


00085-1223 




Hel 


2.1500 


A QQ V 1 ri" 


-15 


-1.04 


34.26 


482.26 






Br7 


2.1661 


l.ZOX iU 


-14 


-2.68 


35.51 


475.85 






HeH 


2.2155 


4.34x10" 


-16 


-0.99 


29.97 


348.24 




Ha 


l-OS(O) 


2.2235 


(S 1 V 1 ri" 


-15 


-1.43 


20.45 


243.83 


00188-0856 


H2 


1-0S(5) 


1.8353 


A 91 V 1 0" 
X lU 


-16 


-1.81 


16.86 


243.82 






Hel 


1.8689 


O.OO X lU 


-16 


-1.79 


18.87 


278.64 






Paa 


1.8756 


5 97x10" 


-15 


-29.4 


27.88 


423.43 




H2 


1-0S(3) 


1.9570 


O.Ul A lU 


-16 


-2.63 


16.45 


215.47 


00582-0258 




Hel 


1.8689 


Q Q/1 V 1 H" 
y .o4 X iU 


-IV 


-1.89 


4.37 








Paa 


1.8756 


1 .00 X LU 


-15 


-31.56 


36.19 


561.35 




Ha 


1-0S(3) 


1.9570 


Z.UDX lU 


-16 


-3.18 


24.26 


294.90 


01077-1707 


H2 


1-0S(2) 


2.0338 


i.zy X iU 


-15 


-1.89 


18.27 


218.47 






Hel 


2.0587 


Z.ZZ X iU 


-15 


-3.27 


26.72 


409.34 




H2 


l-OS(l) 


2.1213 




-15 


-3.45 


23.68 


290.23 






Br7 


2.1661 


5 61x1 0" 


-15 


-9.38 


33.43 


514.02 


02411-I-0354NE 




Hel 


1.8150 


O.oZ X lU 


-15 


-6.95 


23.00 


347.08 




H2 


1-0S(5) 


1.8353 


1.56x10' 


-16 


-2.95 


17.62 


251.22 






Hel 


1.8689 


O.jLiO A _LU 


-16 


-8.58 


26.38 


420.41 






Paa 


1.8756 


O.OU X iu 


-15 


-151.5 


22.67 


339.52 




H2 


1-0S(4) 


1.8920 




-16 


-4.82 


19.11 


305.08 






Br5 


1.9445 


3.29x10" 


-16 


-6.03 


17.39 


272.50 




H2 


1-0S(3) 


1.9570 


3 73x 1 fl" 


-16 


-5.59 


21.31 


301.19 


02411-^0354SW 




Paa 


1.8756 


9 1 S V 1 0" 

Z. To A ±\J 


-15 


-92.3 


24.37 


369.74 


06206-6315 




Paa 


1.8756 


O.UoX lU 


-15 


-30.25 


42.01 


651.92 




H2 


2-lS(6) 


1.8942 


O.DD X iU 


-16 


-3.60 


11.15 


137.46 


19335-3632 




Paa 


1.8756 


U. TtJ A lU 


-15 


-124.0 


24.64 


368.26 






Bt5 


1.9445 


o.iy X iu 


-16 


-7.43 


17.93 


228.89 




H2 


2-lS(3) 


1.9570 


4.07x10" 


-16 


-5.82 


18.79 


264.55 






Hel 


2.0587 


4.71x10" 


-15 


-8.78 


18.23 


237.63 


20100-4156S 


H2 


2-lS(5) 


1.8353 


3.1x10" 


IB 


-11.63 


30.25 


497.00 






Hel 


1.8689 


8.8x 10~ 


16 


-18.61 


62.87 


1008.76 






Paa 


1.8756 


u.uy A lu 


-15 


-136.2 


31.36 


501.56 






Br6 


1.9445 


fi nox 10" 


-16 


-13.91 


35.02 


540.18 




H2 


2-lS(3) 


1.9570 


1.66x10" 


-16 


-39.04 


51.47 


788.84 


20100-4156N 




Hel 


1.8689 


Q oq V 1 ri" 


-16 


-14.08 


33.40 


535.97 






Paa 


1.8756 


1.70x 10" 


-15 


-76.95 


25.46 


407.21 


22206-2715NE 




Paa 


1.8756 


9 RI V 1 n" 


-15 


-37.78 


26.67 


408.40 


22206-2715SW 




Paa 


1.8756 


6 58 X 1 D' 


-16 


-13.18 


28.79 


439.67 


22491-1808 




Paa 


1.8756 


A OA V 1 H" 
4.y4 X iU 


-1!) 


-90.96 


30.07 


449.54 






BrS 


1.9445 


1 91 V 1 0" 
T .ZT X iU 


-16 


-2.00 


12.26 


142.29 




H2 


1-0S(3) 


1.9570 


7 04 V 1 0" 


-16 


-11.93 


34.18 


513.09 




H2 


l-OS(l) 


2.1213 


8 31x1 0" 


-16 


-18.75 


33.77 


459.81 


23128-5919N 




Hel 


2.0587 


0. ( X lU 


-15 


-8.30 


25.75 


427.81 




H2 


1-0S(2) 


2.0338 


O.ZO X lU 


-15 


-18.48 


24.44 


328.09 




H2 


l-OS(l) 


2.1213 


0.00 X lU 


-15 


-13.99 


27.87 


375.23 






Br7 


2.1661 


1.0 / X lU 


-14 


-30.78 


57.11 


811.43 


23230-6926 




Paa 


1.8756 


7.21x10" 


-15 


-111.8 


35.53 


501.97 






BrS 


1.9445 


5.93x10" 


-16 


-8.94 


20.70 


303.25 




H2 


1-0S(3) 


1.9570 


1.30x10" 


-15 


-21.94 


33.92 


497.54 




II2 


i-0S(2) 


2.():!:58 


2.87x10 


l(j 


- i. 15 


i.").i:i 


174.93 


23389-6139 




Paa 


1.8756 


7.34x10" 


-11) 


-77.48 


49.56 


676.14 




H2 


2-lS(6) 


1.8942 


1.71x10" 


-16 


-1.72 


4.65 








Br(5 


1.9445 


9.30x10" 


-17 


-0.90 


2.76 






H2 


1-0S(3) 


1.9570 


2.38x10" 


-16 


-27.34 


43.44 


639.7 






FeH] 


1.9670 


2.63x10" 


-16 


-2.89 


15.24 


150.4 
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uous in the case of a normal star forming galaxy) which, 
in general will diminish the hardness of the ionizing spec- 
trum. In all the following analysis we will make use of the 
recent models for normal star forming galaxies presented 
by Panuzzo et al. I)2UU3|I . In these models the galaxy evo- 
lution is followed with a chemical evolution code and the 
spectro-photometric properties are computed with metal 
dependent stellar population synthesis (Bressan, Chiosi 
&: Fagotto I19H4|I . Line emission is also accounted for by 
means of a large data set of HII region models computed 
with CLOUDY, with the appropriate chemical composi- 
tion. The code (GRASIL, Silva et al. 1998) also account 
for attenuation of light by dust, both in molecular clouds 
and in the more diffuse cirrus component. These models 
have been successfully tested against a number of different 
observations, from UV to the radio (Granato et al. I2UUUI 
Bressan et al. lMT^ and Panuzzo et al. lM75|l . 

We adopt the values quoted by Panuzzo et al. POOS ) 
for a star-forming galaxy, i.e. K^^^'^^ ~3.0, R^"''^"" ~7.2 
and R j„"' ''''~93.5. The latter two values, in particular, 
are lower than that quoted by Hummer & Storey H1987|l 
( ~8.58, and ~103.54), which however refer to a density 
of 10'*cm ~^ and a temperature of Te=10^K. Our sin- 
gle population CLOUDY models show that, already at 
Z=0.02 (solar), the temperature is below 8000K and e.g. 
R °"~8.0. Adopting Rj„"' °" by Hummer & Storey 
H1987|l . would increase the derived E(B-V) by about 0.09 
mag. 

To evaluate the extinction from the observed Paa or 
Br7 fluxes, we have re-scaled our observations to the Ha 
aperture by multiplying our values by the coefficients Cho 
listed in Table 0| As for the extinction law we have used 
that of Calzetti et al. (2000), which is suited for starburst 
galaxies. The derived values of Ay are shown in Table |31 

Table 4. Coefficients for slit corrections. 



IRAS name 


C 


Ha 


CiR 


00085-1223 


1, 


.67 


2.08 


00582-0258 






1.74 


00188-0856 


1, 


.18 


1.44 


01077-1707 


1, 


.28 


3.44 


02411-h0354NE 


1, 


.05 


1.90 


02411-H0354SW 


1, 


.05 


1.88 


06206-6315 


1, 


.32 


2.44 


19335-3632 






3.56 


20100-4156S 


1, 


.21 


1.70 


20100-4156N 






3.35 


22206-2715 


0, 


.41 


3.03 


22206-2715NE 








22206-2715SW 








22491-1808 


0, 


.83 


4.79 


23128-5919 


1, 


.24 


2.04 


23230-6926 


1, 


.16 


3.02 


23389-6139 


1, 


.22 


1.76 



3.3. The star formation rate. 

To derive the star formation rate from the observed line 
intensities we have adopted the calibrations provided by 
Panuzzo ct al. PMjl : 

SFR(Paa) = 5.06x10^41 L(Paa) M0yr-V(erg s'^) 

SFR(Br7) = 6.60x10"'*° L(Br7) M0yr-V(erg s^^) 

For the IR luminosity (Table ^ we have adopted the 
calibration provided by the same authors between the star 
formation rate and the 8/im-lOOO/im infrared luminosity: 

SFR(IR) = 4.63x10-4'' L(IR) MQyr-V(erg s'^) 

All the above calibrations refer to a Salpeter IMF be- 
tween O.IM0 and 120M and take into account an incre- 
ment of 16% due to the lower IMF limit with respect to the 
value of 0.1 5M0 adopted by Panuzzo et al. H2003|l . Paa 
and Br7 fluxes were corrected for extinction and multi- 
plied by the coefficients Cir listed in Tabled Luminosities 
have been computed adopting distances quoted in Table 

m 

The values of SFR, obtained in this way from different 
indicators, are reported in Table IHI The last column of 
Table indicates the ratio between the SFR derived from 
the hydrogen line emission and that derived from the IR 
luminosity. 

On average the Paa luminosity, even corrected for slit 
losses and attenuation, provides a SFR which is only 14% 
of that derived from the far infrared luminosity (assuming 
that the latter is entirely due to the starburst). In the case 
of the objects showing the Br7 line, the average ratio is 
60%. 



4. Discussion 

In all the observed galaxies the Paa or Br7 flux, even 
corrected for aperture effects and extinction derived from 
NIR optical recombination lines (Table EJ, is significantly 
less than that expected from a starburst of corresponding 
bolometric luminosity. A similar "deficit" of recombina- 
tion photons has already been noticed by Goldader et al. 
(1995), based on the analysis of the Br7 line in a sample 
of local ULIRGs. 

Poggianti, Bressan & Franceschini 1)200111 have found 
that the SFR derived from the extinction corrected inten- 
sity of Ha in a sample of ULIRGs was a factor of three 
less than that derived from the FIR luminosity. They at- 
tributed the discrepancy to an age-selective extinction ef- 
fect where a substantial young population was essentially 
escaping the optical detection. On the contrary, in a recent 
analysis of lower infrared luminosity starburst galaxies (10 
<log(L/fl) <11) Mayya et al. l|2004(l have not found strong 
evidence for a loss of optical photons. 

We discuss below several possible explanations of the 
differences between the SFR (and/or attenuation) derived 
from the emission lines and from the IR emission. 
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IRAS name 


Shi3 






Av 










Ha/H/3 


Paa/Ha 


Br7/H/3 


00085-1223 


4.08e-15" 


l.OOe-13" 


6.08 




5.09 


00188-0856 


2.45e-16'' 


3.50e-15'' 


4.52 


4.28 




01077-1707 


1.16e-14" 


9.906-14" 


3.03 




2.97 


02411+0354NE 


8.006-16'' 


4.206-15'' 


1.62 


4.33 




02411+0354SW 


l.lOe-15'' 


7.606-15'' 


2.41 


1.24 




06206-6315 


2.606-16'= 


6.316-15^= 


6.05 


3.25 




20100-4156S 


i.see-is'^ 


1.616-14^= 


3.06 


1.91 




22206-2715 


8.96e-16'' 


6.406-15'' 


2.51 


0.75 




22491-1808 


3.15e-15'' 


2.306-14'' 


2.57 


0.40 




23128-5919 


8.86e-15'= 


8.296-14^= 


3.29 




4.11 


23230-6926 


1.086-15^= 


1.356-14^= 


4.13 


2.39 




23389-6139 


3.40e-16'= 


2.266- 14^= 


8.96 


1.68 





a) Veill6ux et al. |T5^; b) Veilleux, Kim & Sanders itTTOl : c) Due, Mirabel & Maza flMTIl 

Table 6. SFR derived from the observed Paa (or Br7) and IR luminosity. The spectral lines data have been corrected 
for attenuation and for slit losses. 



IRAS name 


L(PaQ) 


SFR(Paa) 


L(Br7) 


SFR(Br7) 


L(IR) 


log 


SFR(IR) 


Line/IR 




[lO^erg s-i] 


[Mq yr-i] 


[lO^'^erg s-i] 


[M0 yr-i] 


[lO-'^erg s-i] 


{Ur/Lq) 


[M© yr-i] 




00085-1223 






4.05 


26.7 


1.08 


11.45 


49.9 


0.54 


00188-0856 


9.1 


48.1 






12.32 


12.51 


570.5 


0.08 


00582-0258 


0.47 


2.4 






2.22 


11.76 


102.7 


0.02 


01077-1707 






8.22 


54.2 


1.92 


11.70 


89.1 


0.61 


02411+0354 


13.22 


66.9 






6.83 


12.25 


316.2 


0.21 


02411-^0354NE 


10.17 


51.5 














02411-F0354SW 


3.04 


15.4 














06206-6315 


8.78 


44.4 






7.73 


12.30 


357.9 


0.12 


19335-3632 


3.20 


16.2 






1.87 


11.69 


86.6 


0.19 


20100-4156 


12.00 


60.7 






20.59 


12.73 


953.1 


0.06 


20100-4156S 


9.08 


45.9 














20100-4156N 


2.92 


14.8 














22206-2715 


6.54 


33.1 






6.49 


12.23 


300.4 


0.11 


22491-1808 


6.38 


32.3 






7.13 


12.27 


330.1 


0.10 


23128-5919 






16.68 


110.14 


3.71 


11.99 


171.8 


0.64 


23230-6926 


15.01 


75.9 






9.29 


12.38 


430.0 


0.18 


23389-6139 


15.42 


78.0 






7.11 


12.27 


329.0 


0.24 



Paa luminosities for IRAS00582-0258, IRAS19335-3632 and IRAS20100-4156N are not extinction corrected, and the SFRs are 
only lower limits. 



4.1. The AGN contribution 

The first and most obvious cause for the photon deficit 
is that a significant contribution to the infrared emis- 
sion comes from an AGN. There is evidence that in the 
most luminous nearby infrared galaxies a significant frac- 
tion of the luminosity could arise in principle from an 
AGN instead of the starburst, as suggested by some opti- 
cal (Veilleux et al. [TM5|l . NIR (Veilleux, Sanders & Kim, 
I1999 t and ISO ( Tran et al. l2001|) spectroscopic studies of 
ULIRGs. 

To explain the observed large discrepancy, the AGN 
should provide a conspicuous fraction of the IR luminosity, 
generally exceeding 80%. 



As already anticipated we have possibly detected 
broad emission line components and [SiVI] emission only 
in IRAS 00188-0856 and IRAS 00582-0258. 

It is interesting to note that, in the IR colour- 
luminosity diagram (Neff & Hutchings [1992), both 
IRAS 00188-0856 with log(S25,™/S60Mm)=-0.63 and 
log(L 6o,.m/L0) = 12.O5, and IRAS 00582-0258 with log(S 
25fj.m/Se0f_im)=- 0.55 and log(L 60Mm/Lo)=11.24, fall in 
the region of overlap between QSO and Seyfert 1 galaxies. 

The presence of a broad component (FWIIM>2000 
km/s) in the Paa line, the possible detection of the high 
excitation [SiVI] line, and their position in the IR colour- 
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luminosity diagram, favour the existence of an obscured 
AGN in IRAS 00188-0856 and IRAS 00582-0258. 



4.1.1. Other diagnostics for unveiling the AGN 

Hints on the nature of the IR luminosity may come also 
from other wavelength observations. 

In star forming galaxies, FIR and radio emissions are 
tightly correlated over a wide range of IR luminosities. 
Sanders & Mirabel ifTMHI report: 



log 



Fpir/{S.75 X IQi^Hz) 
^l.(1.49GHz)/(W m^^Hz" 



2.35 ±0.2 



(1) 



where FFIR=l.26xlO-^\2.58Seo^,„. + ^100^™) Wm'^, 
with S'eo and ^loo expressed in Jy. Radio observations 
of all our target objects have been reported in Tabic ^ 
After extrapolating a few 843MHz data to 1.4GIIz by as- 
suming a radio slope oc i/^"'^, we have evaluated the q 
parameter, reported in Tabic ^ 

AU the objects, including IRASGG188-G856 and 
IRAS00582-0258, but IRAS22491-1808 and IRAS23389- 
6139, fall on top of the FIR-radio relation of starburst 
galaxies. Panuzzo et al. (2003) have shown that in normal 
star-forming galaxies the radio luminosity is an excellent 
indicator of star formation and that deviations from this 
relation arc due to variations of the corresponding FIR 
luminosity as the latter depend more on the details of 
the obscuration. In contrast, in a very young obscured 
starburst there is an excess of FIR emission because core 
collapsed supernovae, thought to provide the 90% of ra- 
dio emission, have a delay of a few million years. For the 
same reason, if the starburst is in a late phase and the 
star formation decreased exponentially, there is an excess 
of radio emission (Bressan, Silva & Granato 12001(1 . 

Based on different arguments. Smith et al. (1998 ) , 
Farrah et al. (i2003p and Fronton et al. (|2004|l . have shown 
that the contribution of the putative AGN to the radio 
luminosity in obscured starbursts (not harbouring radio- 
loud sources) is low. Thus a significant contribution to the 
IR from the AGN would significantly raise the value of q 
above the FIR-radio correlation. This might be particu- 
larly relevant for IRAS00188-0856 and IRAS 00582-0258 
where, as we have already seen, there are hints for the 
presence of the AGN from our own observations, but in 
fact these galaxies show a q parameter typical of star- 
bursts indicating that the contribution of the AGN in the 
IR is not high. 

Another diagnostic for the presence of the AGN is a 
low value of the ratio of the line to continuum emission 
of the 7.7 fim PAH feature, as measured by ISO spec- 
troscopy (Lutz et al. I1998|l . Unfortunately, not all our 
galaxies have this ratio measured. Following this indicator 
(Table three sources could harbour an obscured AGN, 
IRAS 00188-0856, IRAS 23230-6926 and IRAS 23389- 
6139. However only the first source shows evidence for an 
AGN from our NIR observations. This may give further 
support to criticism on the use of this indicator recently 



raised by Farrah et al. (|2003|l (but see also Prouton et al. 
2004). 

There are also XMM-Newton hard X-ray observations 
for IRAS 20100-4156 (Franceschini et al. IMSI- IRAS 
20100-4156 shows a faint X-ray flux of S(2-10 keV)~ 1.9 x 
IQ-^"^ erg s-icm-2, with a (2-10Kev)/(8-1000Afm) flux ra- 
tio ~4.75 X 10~^. Unfortunately, its large distance and 
consequent poor X-ray photon statistics prevented any 
careful X-ray spectral analyses. Although the total source 
flux is slightly above that expected from a starburst of 
similar FIR power (Franceschini et al. I2003|l . there is no 
definite evidence, from hard X-rays, that the galaxy hosts 
an absorbed AGN. For instance MKN231, a typical Syl 
galaxy, is about 34 times brighter, between 2 and 10 keV, 
than IRAS 20100-4156, in spite of having about the same 
IR luminosity. 

In summary, from our own observations and other col- 
lected data in the literature we may conclude that IRAS 
00188-0856 and IRAS 00582-0258 could harbour an AGN. 
However in all cases the contribution of the AGN to the IR 
seems not at the level required to explain the discrepancy 
between line flux and FIR emission. Indeed the first evi- 
dence of a buried AGN in the ultraluminous galaxy UGC 
5101, provided by SPITZER, allows a contribution to the 
total luminosity of only <10% (Armus et al. I2004II . 

4.2. High NIR extinction 

Another alternative to explain the photon recombination 
deficit is simply that the extinction is very high, even at 
NIR wavelengths, and that what we see from optical lines 
is only the external skin of the starburst. 

To evaluate the possible range of that kind of extinc- 
tion we have compared the observed values of the Paa/IR 
(or Br7/IR) ratio with that predicted by models of normal 
star- forming galaxies. In this way we will obtain the ex- 
cess extinction over that of a normal star forming galaxy. 
From Panuzzo et al. ()2003|l we get the following average 
relations between line and IR (8^m- lOOO^m) emission: 

S(Ha)/S(IR)=6. 56x10-3, 

S(Paa)/S(IR)=7.9xl0-4, and 

S(Br7)/S(IR)=6.32xl0-5 

The observed values of S(Paa)/S(IR) or S(Br7)/S(IR) 
are reported in Table [7| The observed intensities of the 
emission lines have been corrected for slit losses by mul- 
tiplying by the factor Cm of Table 01 The average value 
of the observed S(Paa)/S(IR) ratios is ~10% of that pre- 
dicted by the models while, in the case of Br7 line, it is 
~40% of the predicted value. 

Assuming that the difference with respect to the mod- 
els is entirely due to an extinction larger than in the case 
of a normal star forming galaxy, we have derived the corre- 
sponding attenuation in magnitudes, in the line and in the 
visual. The values of Ay must be considered as in excess 
of that of normal star forming galaxies which is typically 
Ay —1 mag. In this way we obtain an average attenuation 
at Paa of about 2.9 mags, corresponding to an average vi- 
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sual attenuation of about 20 mags. These values are much 
higher than those derived from optical emission lines (cfr. 
Table EJ. Values of Ay ~5-50 have been already found 
by Genzel at al. (1998) from mid-infrared spectroscopy of 
a small sample of ULIRGs and are confirmed by recent 
SPITZER observations of selected ULIRGs ( Armus et al. 
[2004,1 . High optical depths (>20) at 1/xm have been also 
inferred from fitting the far infrared to radio SEDs of a 
sample of compact ULIRGs (Prouton et al. I2()()4(l . 



Table 7. Extinction derived from comparison of observed 
and predicted line flux to IR ratios 



IRAS name 


Line/IR 


Obs/Model" 


A(line) 


a'' 






Paa 






00188-0856 


2.73E-05 


3.45E-02 


3.65 


24.9 


00582-0258 


2.13E-05 


2.70E-02 


3.92 


26.7 


02411-f0354 


1.42E-04 


1.80E-01 


1.86 


12.7 


06206-6315 


3.97E-05 


5.03E-02 


3.25 


22.1 


19335-3632 


1.72E-04 


2.18E-01 


1.65 


11.2 


20100-4156 


4.01E-05 


5.08E-02 


3.23 


22.0 


22206-2715 


6.53E-05 


8.27E-02 


2.71 


18.4 


22491-1808 


5.73E-05 


7.25E-02 


2.85 


19.4 


23230-6926 


7.86E-05 


9.95E-02 


2.51 


17.0 


23389-6139 


4.56E-05 


5.78E-02 


3.10 


21.1 












00085-1223 


1.88E-05 


2.98E-01 


1.31 


13.5 


01077-1707 


3.02E-05 


4.77E-01 


0.80 


8.3 


23128-5919 


3.08E-05 


4.88E-01 


0.78 


8.0 



a) Observed L(PaQ)/IR = 7.9E-4 and L(Br7)/IR = 6.3E-5, for 
a normal star forming galaxy (Panuzzo et al. 2003) 

b) AFac«/Ay=0.147 and ABr^/Av=0.097 from Calzetti et al. 
(2000) 



We also notice that, in the case of Paa galaxies, even 
assuming an AGN contribution to the IR of the order of 
50% (which is large, following Prouton et al. 12004(1 . the 
visual attenuation would still remain between 15 and 25 
mags. 

As a second point, we notice that in the galaxies ob- 
served in the Br7 domain the extinction is lower than that 
obtained for the sample with Paa in agreement with a 
slab attenuated model. In fact, for such a model we expect 
A(Br7)~0.66 x A(Paa) (Calzetti et al. 2000). The average 
extinction of the Br7 sample is A(Br7)~0.96, while taht 
of the Paa sample is A(PaQ;)~2.9, with a ratio of 0.33. 
Though we are comparing different galaxies and the statis- 
tics is low (but notice that our Br7 galaxies lie on the re- 
lation defined by the more exhaustive sample of Goldader 
et al. 1997), this may indicate that the photon deficit is 
wavelength dependent^ an thus favour the effects of a com- 
plex extinction geometry over those of absorption by dust 
within HII regions, discussed below. 



4.3. Dust within the HII regions 

In presence of dust within the ionized regions of the star- 
burst, only a fraction / of Ly-continuum photons may 
effectively ionize the gas, while the remaining (1-/) is ab- 
sorbed. This effect causes a deficit of recombination pho- 
tons, with respect to the FIR emission. To explain the 
observed discrepancy between NIR line emission and IR 
fiux, dust within HII regions should absorb about 80% of 
the ionizing fiux. 

Dust absorption within HII regions has been invoked 
by Hirashita et al. 1)2003 11 to explain the anomalous Ha to 
UV flux, observed in a sample of lUE selected star forming 
galaxies. Hirashita et al. (|2,003) have estimated an aver- 
age value of (l-/)~50%, with some objects reaching (1-/) 
~80%. However, it is worth recalling that the conclusions 
reached by Hirashita et al. H2003|l are based on a sin- 
gle screen extinction model. By analyzing the same data, 
Panuzzo et al. H2003|l have instead found that this effect 
may result from "age-selective" extinction, namely that 
younger populations are more extinguished than older 
ones, the latter still contributing to the UV flux but not 
to nebular emission. 

Dust within ionized regions has been indicated also 
by Luhman et al. (2003) as a possibility to explain the 
[CII] deficit relative to FIR, observed in their sample of 
ULIRGs. In fact dust would not only affect the pool of 
ionizing photons, but would also inhibit the penetration 
of 13.6eV— 6eV photons (thought to be responsible of the 
[CII] excitation) in the photo-dissociation regions, while 
preserving the overall FIR emission. 

If confirmed by further studies, a 80% depression of the 
ionizing flux would deeply challenge any determination of 
SFR from even the less extinguished emission lines. 

4.4. Age of the starbursts 

Estimators of the SFR are usually derived assuming a 
continuous star formation rate. This is essentially correct 
when dealing with normal galaxies, but starburst galaxies 
are, by definition, currently dominated by a single episode, 
with a duration which is generally comparable to the char- 
acteristic times of the star formation indicators. In such 
circumstances the use of different calibrators may lead to 
significant discrepancies. For instance, this may be the 
case of the Radio and IR SFR indicators, as discussed by 
Bressan, Silva & Granato H2002|l . 

To get hints on the evolutionary status of our galaxies 
we have analyzed the equivalent widths (EW) of hydrogen 
emission lines (see e.g. Terlevich et al. 2004). These indica- 
tors are relatively unaffected by dust extinction in at least 
two cases, namely when both continuum and line emis- 
sion arise from the same population, or when attenuation 
can be approximated by an intervening absorbing slab. In 
the latter case however, the equivalent width may be af- 
fected by the presence of the old population, outside the 
starburst region. FigurcElcompares the observed Paa and 
Br7 equivalent widths of our galaxies with those predicted 
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Fig. 6. Starburst ages derived from the equivalent width 
of Paa and Br7 emission. The SFR of the models decreases 
exponentially with e-folding time of rsp'fi=25Myr. The 
EW of the models is plotted against the ratio between 
the current age and the SFR e-folding time so that the 
figure does not change significantly by considering other 
plausible values of tsfr- Observed values are plotted as 
horizontal segments. 

by starburst models. The models have been constructed 
by adding nebular emission to the spectral energy distri- 
bution of simple stellar populations of solar metallicity, 
as described in Bressan, Poggianti & Franceschini (j2()()l|l . 
The adopted star formation e-folding time is tsfr='^5i 
but since we are interested in the ratio between the age 
of the starburst and the e-folding time, our conclusions 
do not change significantly by assuming other plausible 
values for tsfr- 

From FigureEl we observe a lack of young objects, with 
only one source, IRAS 23128-5919, appearing younger 
than twice the SFR e-folding time. The other galaxies 
have ages that are either between two and three times 
the SFR e-folding time, or around four times the e-folding 
time. Taken at face value, these estimates would place 
IRAS 00188-0856, IRAS 00582-0258, IRAS 06206-6315 
and IRAS 22206-2715NE in a post-starburst phase. Part 
of the photon deficit discrepancy could then originate sim- 
ply from our use of stationary models. However, given the 
high current IR luminosity of the latter galaxies, their 
bolometric luminosity at the peak SFR would have been 
unreasonably large, between 3x10^'^Lq and 2xlO^*L0. 

It is worth noticing that, being the average distance of 
the Paa sample about four times that of the Br7 sample, 
our slit is sampling an intrinsic galaxy area which increases 
by at least one order of magnitude going from the Br7 
sample to the Paa sample. In spite of that, both the equiv- 



alent width of the lines and the corresponding SFR show 
a comparable distribution in the two samples (though the 
Br7 sample is not statistically significant). This not only 
suggests that the line emission originates from the cen- 
tral regions, but also indicates that the equivalent widths 
are dominated by sources that share the same spatial dis- 
tribution. For instance, they should not be significantly 
affected by an underlying old population, as expected in 
vigorous starbursts (Mayya et al. I2()()4|l . 

We conclude that the observed distribution of EWs, 
skewed toward old ages, is actually another aspect of the 
photon deficit problem, independent from bolometric lu- 
minosity considerations. Unfortunately a statistical argu- 
ment based on the distribution of EWs cannot help to 
discriminate between the different causes. In fact, strong 
age selective extinction could make the starburst to ap- 
pear older than it is; absorption of ionizing photons could 
diminish the equivalent width, because the K band contin- 
uum of the ionizing cluster would suffer much less atten- 
uation; and finally the K-band continuum could be dom- 
inated by the presence of an AGN, still compatible with 
the FIR/Radio correlation (Prouton et al. 2004). 

5. Conclusions and Perspectives 

We have obtained NIR medium dispersion long slit spec- 
troscopy with SOFI at NTT of ten luminous infrared 
galaxies having suitable redshift to push their Paa emis- 
sion in the Ks band. We included also three objects with 
lower redshift, for which the Br7 line was in the Ks band. 

We have found that the Paa emission, even corrected 
for slit losses and for extinction estimated by compar- 
ing our data with optical spectroscopy, is significantly less 
than that expected from a starburst of corresponding bolo- 
metric luminosity. The discrepancy is lower for the galax- 
ies observed in the Br7 region, but in general we confirm 
the existence of a deficit of recombination photons first 
pointed out in ULIRGs by Goldader et al. (1995). 

Furthermore, in IRAS 00188-0856 and IRAS 00582- 
0258 we find evidence for significant broadening of the Paa 
line and for the presence of [Si VI] coronal line. However 
these two galaxies fall on top of the FIR/Radio correla- 
tion, indicating that, though present, the AGN does not 
dominate the far infrared emission. For all other sources 
we do not find evidence of AGN contribution and we ar- 
gue that the studied galaxies appear to be predominantly 
powered by a nuclear starburst. 

Based on current data alone it is impossible to disclose 
the origin of the recombination photon deficit, and we may 
only advance the following hints. 

The galaxies may harbour a highly attenuated star 
forming region. In this case an estimate of the attenua- 
tion can be obtained from the comparison the observed 
linc/IR emission ratio with the predictions of models for 
normal star-forming galaxies. Then the average attenua- 
tion would be Ay ~10-25mag and Apaa — 2-4mag. These 
figures are slightly lower when derived from the Br7/IR ra- 
tio (A Br-y — Imags) and are consistent with the decrement 
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of attenuation going from Paa to the Br7 wavelength re- 
gion. A significant attenuation of the nuclear star forming 
region even in the NIR comes out to be in agreement with 
the large molecular cloud optical depths ( Ti^m > 20), de- 
rived by Prouton et al. I|2U(J4|I from SED fitting of compact 
ULIRGs. It is also in agreement with recent SPITZER 
observations of ultra luminous galaxies with buried AGN 
(Armus et al. 2004). Assuming that the AGN may con- 
tribute about 50% of the IR flux (a quite extreme figure, 
Prouton et al. l2004|l . would imply only a slightly smaller 
attenuation. 

Our finding is also compatible with a scenario where 
a large fraction (~80%) of the ionizing flux is absorbed 
by dust within the HII regions. The required fraction is 
large but comparable with some extreme values found by 
Hirashita et al. (|2003|l in their analysis of a sample of UV 
selected starbursts. This effect has been recently invoked 
as one of the possible causes of the [CII] emission deficit 
relative to FIR observed in some ULIRGs (Luhman et al. 
2003 ). Based on MIR to radio SED fitting, Prouton et al. 
(|2()()4|l have estimated that the dust sublimation radius 
within molecular clouds in compact ULIRGs, is generally 
a fraction of a parsec. Dust may thus survive even in the 
innermost regions of an HII region and, if such a strong 
effect is confirmed, it will challenge our ability to derive 
SFR and/or to study inner environmental conditions, from 
even the most un-extinguished emission lines. 

Though the above alternative scenarios are both com- 
patible with the present data, their predictions at longer 
wavelengths are markedly different and suggest that a de- 
cisive test to disentangle between high nuclear obscuration 
and dust absorption within HII regions, would be to look 
at the Bra 4.05 fim line emission. In fact, assuming the 
Calzetti et al. (2000) extinction law and an intrinsic ratio 
L( Paa)/L(Bra)~3.9 (Panuzzo et al. 12003) we get 

L{Bra)/L{Paa) ~ 0.26 x lo" 04^^. 

Thus if the Paa deficit is due to absorption of ionizing pho- 
tons by dust within HII region and Av—2> mag, as derived 
from optical and NIR-optical emission line ratios, then 
the expected ratio L(Bra)/L(Paa) is ~0.34. On the other 
hand, if the Paa emission comes from a dust enshrouded 
region with say Ay~20 mag, then L(Bra)/L(Paa) ~1.6, 
a factor about five times larger than in the previous case. 
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